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ABSTRACT 

It is proposed that radiation belts similar to the ones in the planetary magnetosphcrc 
can exist for a pulsar with a relatively long period and a strong magnetic field. In the 
belts located in the closed field line region near the light cylinder relativistic pairs are 
trapped and maintained at a density substantially higher than the local Goldreich- 
Julian corotation density. The trapped plasma can be supplied and replenished by 
either direct injection of relativistic pairs from acceleration of externally-supplied par- 
ticles in a dormant outer gap or in situ ionization of the accreted neutral material in 
the trapping region. The radiation belts can be disrupted by waves that are excited 
in the region as the result of plasma instabilities or emitted from the surface due to 
starquakes or stellar oscillations. The disruption can cause an intermittent particle 
precipitation toward the star producing radio bursts. It is suggested that such bursts 
may be seen as rotating radio transients (RRATs). 
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1 INTRODUCTION 

The recent discovery of a new class of radio transient 
sources, known as rot ating radio transients (RRATS) 
(McLaugh lin et aill2006l ). suggests that radio transient phe- 
nomena may be quite common for typical pulsars. About 11 
such sources have been detected so far; they are character- 
ized by short outbursts of a typical from 2 to 30 millisecond 
duration with an average interval between two consective 
bursts from a few minutes to a few hours. Ten RRATs have 
their periodicities determined and three of them have their 
period derivatives identified. Observations appear to sug- 
gest that these objects are part of the normal pulsar pop- 
ulation, with periods of a few seconds, within the period 
range of typical radio pulsars, and magnetic fields up to 
10 10 T, close to a lower limit of the magnetic field of the 
known magnetars. No transient X-rays or optical counter 
part has been observed, though thermal X-ray emission was 
detected from one RRAT, which appears t o be similar to 
X-ray dim isolated neutron stars (XD INS) (Reynold s et al.l 
l200d ; IPopov. Turolla fc Possentill200r3 ). 

Three RRATs with both periods and period derivatives 
identified overlap with normal radio pulsars in the P-P (pul- 
sar period vs period derivative) distribution, which strongly 
suggests that RRATs may have similar properties to radio 
pulsars and should have normal radio emission as well. How- 
ever, regular radio pulses have not been detected from any of 
these sources, suggesting that either (1) these radio bursts 
are distinct from the normal radio emission in aspects of 
emission geometry or processes or (2) normal radio pulses 



are too weak to be detectable l|Weltevrede et aill2006al ). For 
case (1), one could postulate, for example, that radio emis- 
sion jumps between two alternative beaming direction, one 
of which intersects the Earth. For case (2), RRATs may be 
considered as extreme pulses and the 'missing' normal ra- 
dio p ulses can be detected wit h a sufficiently long observing 
time (|Weltevrede et al.ll2006al ). IWeltevrede et al.l (|2006al lbl) 
specifically considered PSR B0656+14, which is a nearby 
pulsar that has pulsed X-rays and intense bursty radio pulses 
similar to RRATs. They suggested that if this pulsar were at 
a distance similar to that RRATs, one would see the strong 
radio bursts as RRATs and its normal radio emission be- 
comes nondetectable. 

If RRATs are indeed pulsars, one has a major diffi- 
culty in interpretation of the observed intense radio bursts 
on the basis of conventional polar cap models because most 
of RR ATs have long periods and pair production is not ef- 
fective. ICordes fc Shannon! |2006) proposed that the bursty 
emission by RRATs is due to circumpulsar asteroids ran- 
domly straying into the magnetosphere; the neutral mate- 
rial is evaporated and ionized, leading to a sudden downward 
flow of charged particles in the polar region, which in turn 
ignites a transient pair cascade toward the PC. The low- 
mass disk hypothesis has some observational support from 
the recent discovery of the fall - back disk of a magnetar, cf. 
IWang. Chakrabartv fc Kaplan! (120061 ) . One main ingredient 
of the Cordes & Shannon's model is that a cascade is as- 
sumed to be due to acceleration in the outer gap. However, 
RRATs are in the regime in the P-P distribution where the 
outer gap is inactive for pair production (through photon- 
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photon collision) except for the case of a nearly orthogonal 
inclination angle where the outer gap may form close to the 
PC and pair production through a single photon decay in 
the magnetic field is important. In this latter case, the cas- 
cade has to be close to the PC and backward emission from 
the cascade may not be able to propagate freely through 
the intervening plasma near the star in the closed field line 
region (CFLR) where induced three-wave interactions are 
expected to be strong (|Luo fe M elrose 2006). 

In this paper, we propose that typical pulsars may 
have radiation belts in the CFLR where relativistic plas- 
mas are trapped. The trapping regions bear many simi- 
larities to the Earth's radiation belts. For example, they 
are similarly subject to various low-frequency disturbances 
that disrupt the trapped plasma causing intense precipi- 
tation. While in the case of the Earth's van Allen belts, 
which are subject to disruption by geomagnetic storms due 
to the solar wind, the proposed main disturbances to the 
pulsar's radiation belts are low-frequency Alfven waves gen- 
erated from the star's surface as a result o f stellar oscillations 
IIMcDermott. van Horn, fc Hansen! 1 19881 ) or shear waves in 
the neutron star's crust due to starquakes (|Blaes et al.l 
1989). It is suggested here that transient radio emission sim- 
ilar to RRATs can be generated as a result of catastrophic 
disruption of the trapping region and that the coherent emis- 
sion is produced by particles precipitating toward the star. 

The density of the trapped plasma can be maintained at 
a much higher value than the local GJ density. In the trap- 
ping region the magnetic field is relatively weak and the syn- 
chrotron decay time is long so that the plasma can be replen- 
ished, leading to a build-up in the density. The major sources 
of the trapped plasma can be particle acceleration in a dor- 
mant outer gap or ablation and ionization of the accreted 
neutral matter. We assume that low-level accretion of neu- 
tral grains from a dust disk or from the ISM occurs and that 
the accreted matter is ionized inside the LC feeding charged 
particles directly to both the open field line region (OFLR) 
llChend[l985l ; lRuderman fc Chendl 19881; ICordes fc Shannon! 
120061 ) and CFLR |Cordes fc Shannonll2006r ). Note that there 
were also previous discussions on accretion of neutral ma- 
terial from the ISM and possible effects on radio emission 
by ioni zation of the neutral matter in the pulsar magneto- 
sphere (lTsvganlll977l ; IWrighdll979l '). Charged particles cre- 
ated from destruction and ionization of neutral material in 
the CFLR can accumulate and be trapped in the region. 
Charged particles created in the OFLR can be accelerated 
in the outer gap and inject relativistic pairs into the trap- 
ping region. High-energy gamma rays emitted by the accel- 
erated particles produce pairs in the trapping region on the 
thermal radiation emitted from the surface. The pair pro- 
duction rate is rather slow, but provided that it exceeds 
the loss rate, it can accumulate plasma near the LC in 
the CFLR. Note that reignition of a lat ent outer gap was 
discussed by iRuderman fc Chend (|l988l ) in the context of 
gamma-ray bursts. However, they considered nearly aligned 
rotators with a period much shorter than that considered 
here. In their model, injection of externally-supplied par- 
ticles in the outer gap is assumed to lead to efficient pair 
creation, a case that does not apply to RRATs or magne- 
tars. 

In Sec. 2, it is argued that plasma trapping regions sim- 
ilar to the Earth's radiation belts can exist in the CFLR of 



typical pulsars with a long period. Injection of charged parti- 
cles to the trapping region due to latent particle acceleration 
in an outer gap or direct ionization of accreted neutral dust 
grains is discussed in Sec. 3 and 4. Sec. 5 discusses the stabil- 
ity of the pulsar radiation belts and possible mechanisms for 
particle precipitation. Application to RRATs is discussed in 
Sec. 6. 



2 RADIATION BELTS 

It is argued here that radiation belts similar to that in plan- 
etary magnetospheres may exist in pulsar magnetospheres 
with relatively weak magnetic fields at the LC. 



2.1 Magnetic mirror 

Particles injected with nonzero perpendicular momenta can 
be trapped in the CFLR due to the magnetic mirror ef- 
fect. For a particle with a perpendicular momentum p±, the 
first adiabatic invariant can be expressed as p\/B = const. 
In this section, we ignore the loss of perpendicular energy 
due to cyclotron decay, which can be important for pulsars 
and is considered in Sec. 2.2. From a quantum mechanical 
view point, this corresponds to a particle remaining in the 
same Landau level. The magnetic field in the CFLR forms 
a 'magnetic bottle' and particles injected at a radial height 
T)i < 1 (rj = t/Rlc is the radius in unit of Rlc = cP/2n) 
can be trapped at rj < 1 provided that their perpendicular 
momenta satisfy the condition: 
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(1) 



where p = (p\ + p 2 ) 1 ^ 2 with py the parallel momentum. 
For example, particles injected with pj_/p|| > 0.35 at the 
LC (r/i = 1) can be trapped; they bounce back and forth be- 
tween two opposite hemispheres in the CLFR above r\ — 0.5. 
The right-hand side can be written into the form sina c = 
(v/Vi) 3 ^ 2 > where the angle a c defines a conal surface in the 
momentum space, referred to as the loss cone at rj. Parti- 
cles with a pitch angle smaller than a c (inside the cone) 
pass through the point rj and those with a pitch angle larger 
than a c (outside the cone) are reflected above rj. For par- 
ticles moving along the last closed field lines, the typical 
bounce time is (cf. Appendix) 



n w -P, 

7T 



(2) 



which is not sensitive to the particle's initial pitch angle on 
at the injection but required to satisfy the condition 
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(3) 



\Rlc ) 

i.e., the reflection radius must be above the stellar radius. 
Eq <[3j) is much smaller than the minimum loss-cone angle 
of the van Allen belts in the Earth's magnetosphere, which 
is a few degrees. However, it is shown in Sec 2.2 that the 
minimum loss-cone angle determmined by cylcotron decay 
can be much larger than ([3]). 

Since the dipole magnetic field is inhomogeneous, 
the guiding centers of trapped particles undergo both a 
gradient drift and curvature drift across the field lines 
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jNorthrop fe TelleiT l960). These drifts cause particles to cir- 
culate around the star forming a ring current, similar to the 
ring current in the Earth's radiation belts. The current gives 
rise to perturbations to the dipole field at 7] = 1: 



«l(T) 
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(4) 



SB NLmeC 1 ^ 
~B~ L ~ 2Ubl 

where Ubl = B^/fio is the magnetic energy density at the 



1/2 



LC, N L = JV GJ 6I2, d = {Ro/RlcY* 1 « 0.012(P/1.5s 
is the half-opening angle of the PC, and Ngj ~ 2.3 x 
10 17 (B s /10 9 T)(P/1.5 s)" 1 m -3 is the GJ density at the PC 
where the magnetic field is assumed to be B 3 . In practi- 
cal situations where the plasma density can be much higher 
than the GJ density (Sec 3), the perturbations can exceed 
the estimate given by Since the ring current can vary 
temporally as a result of instabilities in the trapped plasma 
(cf. Sec 5), SB should be time dependent. The effect of such 
perturbations on pulsar electrodynamics and their possible 
contribution to pulsar timing noise will be discussed else- 
where. 



2.2 The loss-cone angle 

The lifetime of trapped particles is limited by cyclotron de- 
cay and pitch-angle scattering, both of which violate the first 
invariant. The electron's synchrotron cooling time, i.e., the 
time that it loses half of its energy (m e c 2 7), is t s = T0/7 
with 

(5) 



TO = 



where r e w 2.8 x 10 -15 m is the classical radius of the elec- 
tron and fi e is the electron cyclotron frequency. Figure [1] 
shows a pulsar distribution on P-P with the time ([5]) indi- 
cated. Six examples where t is calculated are listed in Table 
1. For fast rotating young pulsars or millisecond pulsars, ((5]) 
is very short ~ 10 -6 s, and for typical pulsars it can be quite 
long~ 10 4 -10 5 s. Since in the synchrotron regime, the pitch 
angle remains constant, the relevant time that constrains 
particle trapping is the cooling ti me in the limit p± —> 1, 
7 -> 7c = v%[l+pi|o/(l+pio)] 1/2 (|Luo. Shibata fc Melrose! 
2000), where p±o and p« are the perpendicular and parallel 
components of the particle's dimensionless initial momen- 
tum. This gives t s ~ to/^c, which has the following numer- 
ical form 



(6) 



One has t 3 ~ to if a particle is injected into the trapping 
region with a large initial pitch angle tanao = P±o/P||o ~ 1- 
However, one has a much shorter cooling time if the initial 
pitch angle is small ao <C 1. In the cyclotron regime, the par- 
ticle radiates at a much slower rate, relaxing to the ground 
Landau state on a time T a jc ~ toj c - 

Equating the cyclotron decay time to the bounce time, 
one finds the minimum pitch angle for a particle to be re- 
flected at the mirror point. This angle is estimated to be 

1/2 / p x -1/2 



on 5s 0.04 



B s 



10 9 T 



(7) 



For particles to remain trapped for time much longer than 
rs, the pitch angle at injection needs to be much larger than 
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10 s 



P(s) 

Figure 1. Pulsar distribution in P-Bl with the decay time 
to shown. The magnetic field at the LC is defined as 
= B a 0®. The pulsar sample is from ATNF pulsar cat- 
alogue (http://www.atnf.csiro.au/research/pulsar/psrcat). The 
three large shaded dots represent RRATs. The lighter square rep- 
resents a typical magnetar. 



Pitch angle scattering by waves can lead to particle dif- 
fusion to small pitch angles and limit the life time of the 
trapped particles. Since the decay time decreases rapidly 
with decreasing altitude according to r s oc rf , particles scat- 
tered to small pitch angles can reach low altitudes rapidly 
radiating away their perpendicular energy and eventually 
reach the star. For example, particles with f c ~ 1, which 
are scattered into the loss cone of ct c = 0.01, can travel 
down to an altitude r\ = 0.05^ where the decay time is 
only t s ~ (lOOms)??;. The sources of the waves can be inter- 
nal, i.e., they are generated in situ through plasma instabil- 
ities, or external, i.e., they are generated outside the region. 
Since particles with small pitch angles escape through both 
ends of the 'magnetic bottle', there is an excess of particles 
with a large pitch angle; this gives rise to an ideal condi- 
tion for instabilities to develop (driven by an inversion in 
the particle pitch angle distribution). Waves generated from 
instabilities can cause pitch angle diffusion, further enhanc- 
ing particle precipitation. In principle, this process can im- 
pose an upper limit to th e density of the trapped plasma 
IlKennel fc Petschecldll966h . However, since the efficiency of 
pitch-angle scattering depends explicitly on the wave inten- 
sity that in this case is limited to much smaller than the 
density of the particle kinetic energy, one needs a relatively 
high plasma density for such upper limit to be significant 
(cf. Sec 5). Alternatively, pitch-angle scattering can be sig- 
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nificant if waves that propagate into the region can interact 
with the trapped particles in resonance. This possibility is 
discussed in Sec. 5. 



2.3 Sources for the trapped plasma 

Particle loading into the trapping region can be achieved 
through a process that is not constrained by the first adi- 
abatic invariant. We consider two possible sources for the 
trapped plasma, both of which involve extrinsic charged par- 
ticles: (1) pair production in the trapping region by high- 
energy gamma rays emitted by particles accelerated in an 
outer gap and (2) direct ionization of neutral material in 
the thermal radiation field in the CFLR. Both cases may 
involve ionization of neutral matter migr ating into the mag- 
netos phe re du e to accretion of a disk (|Cordes fc Shannon 
20061 ; 120061 ) or neutral dust grains in the ISM (|Chene 
198a ). For (1), the outer gap, which is usually dormant for 
typical pulsars, is reignited with the supply of charged par- 
ticles. Although the gap is rather inefficient for pair produc- 
tion, it can lead to a build-up of particles in the CFLR over 
a time much longer than the pulsar period (cf. Sec. 4). For 
(2), charged particles from the ionization can be accelerated 
through either inward (toward the star) cross-field diffusion 
or resonant wave-particle interactions, which is similar to 
the particle acce leration in the van Allen belts in the Earth's 
magnetosphere (Ho me et al.l 12005). 

In principle, high-energy cosmic rays can also inject 
pairs into the trapping region. (Note that cosmic rays con- 
tribute to trapped plasmas in the van Allen belts through 
decay of upward deflected neutrons produced by high-energy 
cosmic rays.) For example, cosmic protons with a gyroradius 
larger than Rlc can drift into the pulsar magnetosphere. 
For a typical pulsar, the proton needs to have a Lorentz 
factor 7 > Qp/Q. ~ 5 x 10 4 , where Q p is the proton cy- 
clotron frequency at the LC and £1 = 2n/P. Ultrarelativis- 
tic protons produce pairs on the thermal radiation through 
p + 7 — > p + e . However, it can be shown that such mecha- 
nism is not effective. Assuming a typical proton flux density 
similar to that received on the Earth, F p ~ 5 x 10 -7 m -2 s _1 
at 10 5 GeV, one finds a pair injection rate about 10 6 s _1 . It 
would take much longer than the Hubble time to fill the 
magnetosphere to the GJ density. 
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3 DORMANT OUTER GAP 

We show that a dormant outer gap can supply pairs to the 
trapping regions. It is generally thought an outer gap exists 
between the null surface and th e boundary surface sepa- 
rating the OFLR and the CFLR (ICheng. Ho. fc Rudermanl 
1986 1; iRomanil 1 19961: Iffirotani fc Shibatall2002l ; IZhang et"aLl 
2004 ; iTakata et al.ll2006h . For young, fast rotating pulsars, 
particle acceleration in the outer gap can self-sustain a pair 
cascade, which can transfer a significant fraction of the gap 
potential into the secondary particles kinetic energy, which 
is then radiated at high energies. As the pulsar period in- 
creases, the gap become less efficient in pair production and 
thus, for slow pulsars, pair production in the outer gap be- 
comes negligible. However, provided that there is an exter- 
nal source of particles, the potential energy of such dormant 




Figure 2. Radiation belts (shaded area) and transient radio emis- 
sion (TRE). Radio bursts are produced due to particle precipita- 
tion toward the star in the closed field line zone. The thick line ab 
along the last open field line is the outer gap where electrons are 
accelerated toward the star. The curvature photon (7) emitted 
by these electrons can be converted to a pair by colliding with 
thermal photons emitted from the star's surface. Pairs are then 
trapped in the shaded region due to the magnetic mirror effect. 
Disturbances to the region, such as Alfven waves launched from 
the star, can lead to an intense plasma precipitation. As particles 
move toward the star, coherent radio waves can be generated, 
e.g., through synchrotron/cyclotron maser or two-streaming in- 
stability. 



gap can be tapped and transferred to plasma kinetic energy 
trapped in the CFLR. 



3.1 Acceleration in the outer gap 

The accelerating electric field along the field lines in an outer 
gap can be written as a fraction of the maximum potential 
drop <E> m across the PC: 



Rlc 



(8) 



where w g ^ 1 is a parameter characterizing the ge- 



ometry of the gap and $„ 



0.56%B s Ro 



1.7 x 



10" V(B s /10 a T)(P/2s) . Since pair production can 
change the gap geometry, the parameter w g is usually de- 
termined self-consistently by taking into acco unt the effect 
on th e accelerating potential by pair creation |Takata et alj 
2006). As we consider here specifically slowly rotating pul- 
sars, in contrast to rapidly rotating young pulsars such as 
the Crab pulsar, pair production has little effect on the gap 
electrodynamics and thus, one can treat w g as a constant. 

Assume that the outer gap is located along the last 
open field line starting from the null surface at a radius rj NS 
and consider acceleration of an electron toward the star. 
Possible extrinsic sources of charged particles are discussed 
in Sec. 3.3. When curvature radiation is the dominant en- 
ergy los s process, an assumption used in the usual oute r gap 
models (|Cheng. Ho. fc Ruderman1ll986l ; lRomani|[l996T ). the 
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radiation-reaction limited Lorentz factor is 

'Pi C ^ 1/4 



7c 



= 2 



1/2 9 7™??ns) 1/4 (~^f~) 

1n 7 1/4 1/4 / P Y 1/4 ( B 3 X 1 / 4 



4.1 x 



where 7™ = e$ m /m e c 2 « 3.2 x 10 7 (P 3 /10 9 T)(P/2 s)~ 2 , 
Vns = 0.5?7 NSl o.5, the curvature radius at ri NS is assumed 
to be R c = (4/3)77^^0(1 - 3t7 ns /4) 3 / 2 /(1 - % s /2), and 
r e ~ 2.8 x 10 15 m is the classical radius of the electron. The 
characteristic energy (in m e c 2 ) of curvature photons emitted 
by these electrons is 



10 9 T 



\3/4 

) ,(10) 

where 7 = min{7 c , w g ^/ m }, the approximation is obtained 
with 7 = 7c, and A c = Ti/m e c w 3.9 x 10 -13 m is the Comp- 
ton wavelength divided by 27r. One may estimate the pulsar 
period range where the acceleration is limited by curvature 
radiation: 

3/7 

(11) 



-1/7 f 



P< 1.8w 4/7 »7 NSi0 . 5 Uo9T 

The probability of a photon being converted to a pair 
through collision with thermal photons is tiny, and hence 
the bulk of curvature photons escape as a MeV gamma- 
ray flux with a luminosity L 1 « e c NLCTv(r] ¥SS RLc) 2 'm e c 2 ~ 
10 19 Js _1 . For slowly rotating pulsars this flux may be too 
low to be detectable. 

For pulsars with rotation periods that do not satisfy 
(llip . primary particles may lose energy principally through 
inverse Compton scattering (ICS) on thermal radiation from 
the star's surface. For a surface temperature T s , the number 
density of thermal photons with a characteristic energy }cbT s 
is 

3 / D N 2 



\vRlc 

T, 



7T 2 V A J 



(12) 



where T = k B T s /m e c 2 « 1.7 x 10~ 4 (T S /10 6 K). Since the 
thermal photons propagate radially, backward moving elec- 
trons undergo head-on scattering on these photons. In the 
electron rest frame, the thermal photon energy is 7OT 3> 1, 
implying that the scattering is in the Klein-Nishina (KN) 
regime. The energy loss rate is 



d 7 ca T 0jOl 



In (2 7 e 2 



(13) 



dt 16A! n 2 
with ot the Thomson cross section. 

3.2 Pair creation 

The only viable channel for pair creation far from the star's 
surface is through photon-photon collision. The primary 
photons can be produced from curvature radiation or ICS. 
First, consider pair production due to curvature radiation. 
The number of curvature photons with energy e > e c emit- 
ted by a primary electron over a distance AqRlc is written 
as 



3 3/2 a/ Ag7 



1/2 

??NS 



(14) 



I(e/e 



dy / dtK 5/3 (t), 



(15) 



where Ag < 1 is the thickness (in units of Rlc) of the 
shell region, otf m 1/137 is the fine constant and K n (t) is 
the modified Bessel function. The number of curvature pho- 
tons emitted decreases rapidly with increasing e/e c . One has 
n c « 6.4 x 10 3 for e ^ e c and n c « 1.7 x 10 3 for e > 2e c 
and n c w 2 x 10 2 for e > 4e c , with Ag = 0.2, 77 NS = 0.5, and 
7 = 3.2 x 10 7 . The number of pairs per electron produced 
near the LC, i.e., the multiplicity denoted by Ml, is esti- 
mated as Ml ~ ncUphajjRLcAg where er 77 w 2 x 10 -29 m 2 
is the cross section of pair production via photon-photon 
collision. Since pairs are produced on the Wien tail of the 
thermal spectrum, in stead of (|12|l , the thermal photon num- 
ber density is written as n p n ~ (2/7t 2 j7 2 )(0t/A c ) 3 £w, where 
€,w = (eth/Or) 2 exp(— eh/9t) is the parameter characteriz- 
ing the Wien tail (e p h 3> Ot) where thermal photons satisfy 
the pair production threshold, 



(1 - cos<9)e 



(16) 



where 9 is the angle between the pair-producing gamma- 
ray and thermal photon. One has £w ~ 0.06 for e — 4e c , 
6 — tv / 3 and the numbers given in QD) . Assuming that |TT} 
is satisfied, one obtains 



4.3 x 10~ 3 w 



V0.06/ 



1/4 ( t s \ 3 (_p_y 5/4 ( b s y/ 4 

9 U0 6 K/ V 1.5s J U0 9 T/ 



Ago. 2 



(17) 

with Ago. 2 = 0.2Ag. The multiplicity increases rapidly with 
decreasing pulsar period. As examples, from Eq (|17[) one 
obtains M L « 0.027 for J1913+1333 with T s = 10 6 K and 
Mi w 0.015 for J1819-1458 with T s = 1.4 x 10 6 K, where 
i) NS = 0.5 and Aq = 0.2 is assumed. 

Pair production due to ICS becomes significant if the 
curvature photon energy is too low to satisfy the threshold 
condition (|16[) . Since the scattered photon energy is e 7 ~ 7 
in the KN regime, the production rate of the scattered pho- 
tons with energy e 7 is dn ics /dt as —y~ 1 d'y/dt. The number 
of photons emitted through ICS over a distance ?7 NS Pic is 



gtRlc g|g| ln(2 7m O T ) 
16A 3 rj NS 7 m 



(18) 



Similarly, one has Mi m n IC 
the following numerical form 



^Upha-y-yRLcAg, which takes 



M L 



2.1 x 10" 



uo 6 K7 



Bs 



(io" 1 



Ago. 2 
n 3 

'/ NS ,0.5 



(19) 



Pair production is strongly dependent on the surface tem- 
perature but insensitive to the pulsar period. 

In the outer gap models, the primary particles are 
assumed to be created in a pair cascade in the gap 
ijCheng. Ho. fc Ruderm"anlll986l : lRomani|[l99r3 l and the pri- 
mary flux is limited to a value near the GJ flux. To allow a 
substantial energy transfer from primaries to pairs, as this 
is the case for high-energy pulsars, one must have the pair 
multiplicity in the vicinity of the gap to be at least one, i.e., 
M g 1, with M g = 1 defining a deathline where the gap 
becomes inactive. For parameters appropriate for RRATs 
and magnetars, one has M g <C 1 except for a — > 7r/2. Thus, 
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Table 1. The magnetic field Bl (at the LC), Nl = Nqj9®, synchrotron decay time T d times 7 
and charge loading time tl- The magne tic field on the polar cap s urface is obtained from B s = 
6.4 x 10 15 (PP) 1 / 2 l IUsov fe Melroselll995h . twice the value quoted bv lMcLauehlin et al.l ll2006ri . The 
time tq is calculated at 0.5Rlc- 



RRATS/Pulsars 


P (sec) 


B s (10 9 


T) 


B L (10- 4 T) 


N L (m" 3 ) 


TO 


(sec) 


t l (sec) 
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one can ignore the feedback effect by pair production on the 
accelerating electric field in the gap. 



3.3 Source of primary particles 

Low-level accretion of neutral grains from a dust disk or the 
ISM can provide extrinsic charged particles. The accreted 
neutral material can be ionized inside the pulsar magneto- 
sphere by the thermal radiation from the star's surface and 
particles produced from ionization are then channeled to the 
outer gap or trapped in the CFLR. The maximum flux that 
can be tolerated in the gap is about the GJ flux cNl/i]^ s 
with the gap location assumed to be along the last open field 
lines starting from the null surface at the radius »7 NS . It is 
unlikely that the supply of charged particles matches exactly 
the GJ flux and the acceleration may become oscillatory so 
that the net flux is maintained at about the GJ level. For the 
case of accretion of the ISM dust grains, the accretion rate 
is M ~ {GM) 2 p IS M/v 3 ~ 10 4 kgs~\ where G is the gravi- 
tation constant, M is the neutron star's mass, and v is the 
pulsar velocity, Pism is the ISM density. Dust grains of small 
size (< 0.1/im) may not survive sputtering by protons when 
crossing the bow shock and those with size ~ 0.1/4 m have 
a better chance to reach the LC (C heng|[l985h . Assuming 
that flux rate F e of electrons into the gap is 10% of this flux 
rate, one estimates F e = cN a ~ 0.1(c/2ttR lc )(M /m p ) » 
5 x 10 13 m _2 s _1 . Here we assume v = 100 kms -1 and 
Pism = 10~ 21 kgm -3 . One finds this flux is considerably 
larger than the GJ flux cN L ~ 2.4 x 10 13 m -2 s _1 . The ac- 
celerati ng electric field changes si gn if the flux exceeds the 
GJ flux (jHirotani fe Shib ata 2002). Although a div ergent so- 
lution for E\\ is possible ( Hirotani fe Shibatall2002T ). here we 
assume that the maximum accelerating potential is limited 
by ([8j determined by the vacuum potential drop. 

For pulsars with T B > 10 K, the ISM grains may not 
survive disintegration |Cheng||l985l ) and in this case the disk 
is the main source of the accrete d neutral matter. The dis k 
model was discussed in detail bv lCordes fe Shannon! |2006) . 
The recent discovery of a fall-back disk around a magnetar 
lends support for the existence of such disks around long- 
period pulsars. 



3.4 Particle loading in the CFLR 

Despite the slow rate of pair production, a significant num- 
ber of particles can be accumulated in the CFLR over a time 



2> P. A schematic of pair injection by a latent outer gap is 
shown in Figure [5] 

One may estimate the time required for the plasma den- 
sity to increase to the GJ density in the trapping region. 
Assume a particle flux in the gap is cNl / ??^ s . Acceleration 
of these particles can inject pairs into the CFLR at a rate 
itMNlcS s Rlc/Vns w hh the gap's cross section area approx- 
imated by tSsTi^Rlc (where 5 S < 1). A continuous injec- 
tion of pairs can lead to build-up of the plasma density in 
the trapping region. The time for the density to reach the 
GJ density Nl is 



TL 



AgP 

2ttMS s 



(20) 



We refer to (|20|) as the pair loading time. Numerical exam- 
ples of tl are shown in table 1, where we assume A^i = 0.2, 
i) NS = 0.5, and 9 = 7r/3. The surface temperature is as- 
sumed to be T s = 10 6 K for J1913+1 333, T s = 1.4 x 10 B K 
for J1819-1458 (|Revnolds et al.ll2006l ) and T 3 = 1.5 x 10 B K 
for a magnetar. One has r c ^> tl for these three pulsars. 
For J1317-5759, one finds t l ~ 1.5 x 10 5 s for 9 = tt/2, 
T a = 1.5 x 10 6 K and <r) ns = 0.2. 

Alternatively, the plasma can be injected near the star 
throu gh photon decay in the magne tic field |Wang et al 



1998). Such case was also considered bv lRafikov fe Goldrcich 



(2005) for the double pulsar system. In their model, they 
considered pitch angle increase due to absorption of the co- 
herent radio waves from the other pulsar in the system. 



4 INJECTION THROUGH IONIZATION 

Radiation belts can form directly through evaporation and 
ionization of neutral grains accreted from a disk or the ISM. 
One can show that the density of trapped plasmas can eas- 
ily exceed the local GJ density with a slow accretion rate. 
For example, to accumulate a plasma with the GJ density 
over a time tl, a rate as low as cNl/tl would be sufficient. 
The trapped particles can be accelerated to relativistic en- 
ergy through either cross field diffusion toward the star or 
interactions with waves near the cyclotron resonance. Both 
processes are believed to be responsible for particle accel- 
eration in the Earth's van Allen belts where electr ons are 
accelerated up to 10 MeV energy (|Horne et al.ll2005l ). In the 
cross field diffusion, extremely low-frequency disturbances 
(ELFWs) in resonance with the particle drift frequency 
(around the star) can force particles to diffusion across the 
field lines inward (toward the star) where the magnetic field 
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increases. Conservation of both the first adiabatic invariant 
pj_/B and the second adiabatic invariant <f> dsp\\ (bouncing 
back and forth between two mirror points) implies an in- 
crease in both p± oc 1/r 3 and pii oc 1/r; thus, p increases 
as tana = p±/p\\ oc 1/r 1 ^ 2 increases with decreasing r (in- 
ward), that is, the particles are accelerated. For the Earth's 
radiation belts, ELFWs are generated due to the large scale 
spatial disturbances by the solar wind. In the pulsar case, 
the most plausible candidate for such low frequency waves 
is Alfven waves (cf. Sec. 6). 



5 STABILITY OF THE TRAPPED PLASMA 

The pulsar radiation belts have many similarities to plane- 
tary radiation belts. For example, in both cases the stability 
of the trapped plasma can be strongly affected by change 
of the global magnetic field structure, i.e. distortion of the 
'magnetic bottle', or strong pitch-angle scattering of parti- 
cles on waves. For the Earth's van Allen belts, such activi- 
ties are due to the disturbances from the solar wind, while 
activities in the pulsar radiation belts are powered by the 
rotational energy or magnetic energy. In magnetars, strong 
starquakes can lead to significant distortion of magnetic field 
lines, but such major events may not be common for nor- 
mal pulsars. Here we only discuss low-level starquakes or 
neutron star oscillations (cf. Fig [3}, which can lead to pitch- 
angle scattering by Alfven waves generated at the surface. 



5.1 Pitch-angle scattering 

A formal theory for diffusion of particles in momentum space 
is the quasilinear diffusion formalism in which the processes 
arising from a feedback from growth of linear waves are 
determined by diffusion coefficients. In momentum space, 
such diffusion can be separated into diffusion in pitch an- 
gle and in p. The general formalism is outlined in Appendix 
B. Here we consider specifically the pitch angle diffusion 
coefficient D aa (cf. Appendix). Quasilinear diffusion occurs 
effectively through resonant wave-particle interactions. The 
most widely-discussed process in the context of cosmic ray 
confinement is the scattering due to c yclotron reson ance 
w - fe||W|| - sfi e /7 = with s = ±1 <|Melrosdll980l ). In 
this case, the wave fluctuations can be regarded as predom- 
inantly magnetic and as the zeroth order approximation in 
the expansion on uj /kc <C 1, the process can be viewed as 
elastic scattering of particles by magnetic fluctuations, caus- 
ing change of their motion direction. Since the freque ncy 
of waves generate d from starquakes dBlaes et al.l Il989l) o r 
stellar oscillations l|McDermott. van Horn, fc Hansen Il988l ) 
is low, about kHz as compared to the cyclotron frequency 
> 10 MHz, the cyclotron resonance condition is generally not 
satisfied except for the cases where these waves can cascade 
into high-frequency waves, e.g., through three- wave interac- 
tions. Here we only consider pitch-angle scattering due to 
the Cerenkov resonance. 

Unlike wave-particle interactions in the cyclotron reso- 
nance, the Cerenkov resonance only changes p« not p±. In 
general, it leads to an increase in pu , corresponding to a de- 
creasing pitch angle a. There are two relevant low-frequency 
wave modes: the X mode and the LO, corresponding re- 
spectively to the fast mode and the Alfven mode in the 



MHD. Pitch angle diffusion due to th e fast mode waves in 
Ceren kov resonance was discussed by ISchlickeiser fc Miller! 
(1998). Since there is a nonzero SB along the mean mag- 
netic field, fast particles in the Cerenkov resonance bounce 
back and forth between two successive magnetic compres- 
sions (acting as two magnetic mirrors) in the wave frame. 
On average, particles gain energy causing an increase in pu . 
This process, also called transit-time damping, works for a 
low Alfven speed va/v <C 1, a condition not satisfied for 
pulsars. So, instead we consider a low frequency LO mode, 
which we refer to as the Alfven wave. In kinetic theory in 
the limit va — > c, which is appropriate for pulsar magneto- 
pheric plasmas, Alfven waves have a nonzero SE along the 
magn etic field (|Arons fc Barnardl Il986l ; iMelrose fc Gedalinl 
Il999t ) and particles can be accelerated through the Cerenko v 
resonance (|Volokitin. Krasnosel'skikh. fc Machabelil Il985l ). 
The scattering time t a = 1/ D aa may be written in terms of 
the ratio of the magnetic energy density Ub = SB\/2hq of 
the wave to the plasma kinetic energy density U p = m e c 2 (7) 
(Appendix) 



1 / Up \ /Up 
2trj V U B J \w 



sin a cos a 



(21) 



Considering a plasma consisting of electrons and positrons, 
for Ub ~ Up, u = O.luip, 6 — 0.3 and a = 7r/4, one estimates 
t a ~ Is. One concludes that diffusion due to particles in 
Cerenkov resonance with a low frequency Alfven wave can 
effectively transfer particles to small pitch angles. 

It is appropriate to comment here that growth or damp- 
ing of Alfven waves through the Cerenkov resonance is gen- 
erally not effective because the parallel component of the 
wave polarization is rather small, en ~ (u/uip) 2 <C 1. Thus, 
quasilinear diffusion as a result of such linear wave growth is 
not significant, limited by Ub /U p <C 1 (cf. Eg. I21[) . This has 
been the main reason that such two-stage diffusion mech- 
anism is not favored for the interpretation of cosmic-ray 
propagation, which, as inferred from observations , is sub- 
ject to strong-pitch angle scattering in the ISM (| Melrose! 
ll980T ). However, such limit is not applicable in our case as 
Alfven waves are assumed to be generated externally, not 
through the Cerenkov resonance, with magnetic energy den- 
sity comparable with U p ; thus, efficient diffusion occurs even 
for {u/iOp) 2 < 1. 



5.2 Precipitation 

Disturbances to the radiation belt can be catastrophic if the 
wave that propagates to the region is in the form of short 
bursts producing Ub ~ U p in the trapping region. As shown 
from (|21[l. the scattering time can be as short as seconds. 
The requirement of Ub ~ U p is rather modest. For example, 
for (7) = 5, one has U p ~ 5 x 10~ 7 J m -3 , giving a luminosity 
AgR 3 LC U p /At w 10 16 Ago. i(l s/At) W where At is the pulse 
duration. As this luminosity is minuscule compared with 
transient emission in magnetars (typically 10 25 W), these 
waves are not significant in producing high-energy emission. 
Strong pitch angle diffusion channels most particles into the 
much smaller loss cone through which particles can reach a 
region close to the star where their synchrotron decay time is 
short. This can cause a sudden intense precipitation of parti- 
cles on a time scale shorter than the pulsar period. One may 
estimate the required wave amplitude in terms of magnetic 
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Accretion+ionization 



pairs+ions \ Radiation Belt 



GJ background plasma 



Aflven waves (A) 




7 from an outer gap 



Figure 3. Disruption of the trapping region and particle precip- 
itation in the CFLR. Alfven waves (A) emitted from the surface 
can bounce back and forth in the closed magnetic flux tube in- 
teracting with the trapped particles. LS is the direction in which 
an observer would see the bursts arising from streaming instabil- 
ity as a result of downward moving particles. In contrast, bursts 
arising from cyclotron maser at reflection points can only be seen 
in the direction approximately perpendicular to the field lines. 



fluctuations as SBa/Bl ~ (U p /Ubl) 1/2 ~2x 10" 4 , imply- 
ing low intensity quakes or oscillations would be sufficient 
for triggering disruption of the trapped plasma. 

If the triggering mechanism is starquakes, the reoccur- 
rence time of such particle storms is essentially the reoc- 
currence time of the quakes. The physics of neutron star 
quakes or oscillations is not well understood. A possibility 
is that a strong toroidal field exists in the crust underneath 
the surface and its stress can be released as quakes trans- 
fering the magnetic energy to the magnetosphere in A lfven 
waves (IBlaes et al ] 119891 ; iThompson fc Dunxjj Ejjl). One 
expects that such starquakes are relatively more frequent 
for a high-field, young pulsar. If generation of low-intensity 
Alfven waves is so frequent that the reoccurrence time is 
much shorter than the pulsar period, trapped plasmas do 
not have enough time to accumulate and radiation belts may 
not form. This may be the case for magnetars with active 
transient emission. 



6 APPLICATION TO RRATS 

An important consequence of the formation of pulsar radi- 
ation belts is transient emission — radio bursts that can be 
produced as particles precipitate and such transient emis- 
sion may be relevant to RRATs. 



6.1 Coherent emission 

RRATs are bright radio bursts with a flux density up to 
a few Jy at 1.4 GHz. As for normal radio pulses, the 
emission mechanism must be coherent. Since there is no 
widely accepted mechanism for the pulsar radio emission, 
it is expected that many uncertainties would remain in 



identifying the specific emission mechanism for the tran- 
sient case. Nonetheless, in the following we consider two 
possibilities in particular: two-stream instability and cy- 
clotron/synchrotron maser. 

Two-stream instability 

The two-stream instability has been discussed extensively 
as a possible mechanism for pulsar radio emission, but the 
growth rate is not sufficient for the instability to develop in 
the polar region, mainly because o f the very h igh Lorentz 
factor of the primary particle beam (|Usovlll987r ). It is shown 
here that the plasma conditions created during the particle 
precipitation is favorable for such an instability to occur. 

Two-stream instability can develop when particles tra- 
verse a corotating background plasma which is assumed to 
be stationary. Assuming the precipitating plasma density is 
2Ml times the local GJ density Nlii~ 3 and the spread in 
the particle momentum distribution is ignored (this spread 
can reduce the growth rate, cf discussion below). Consider 
wave growth due to the Cerenkov resonance u) « k\\V\\. The 
characteristic frequency is determined by the denser plasma. 
For 2Ml <C 1, the characteristic frequency is the frequency 
uip = (e 2 NL/eom e ) lll2 rj~ 3 ^ 2 of the stationary background, 
corotating plasma, with the backflowing plasma regarded as 
a weak beam satisfying the resonance condition. The growth 
rate is then estimated to be 



2 



M 



1/3 



2-7 



(22) 



The growth rate reaches maximum at 7 — > j c in the limit 
p± — > 1. For 2Ml 1, the growth rate can be estimated 
in the rest frame of the backflowing plasma assuming that 
the corotating background plasma is a beam in Cerenkov 
resonance. One finds 



(4M L )V3 7 ' 



(23) 



where u> = u! p (2Mi J 'y) 1 ^ 2 is the characteristic frequency seen 
in the pulsar frame. 

Note that for two-stream instability in the polar cap 
model, the energetic particle beam is usually identified as 
those primary particles from the PC, with a Lorentz factor 
as high as ~ 10 6 — 10 7 and i nstability is strongly suppressed 
l|Cheng fc Ruderman|[l977l ; IUsovlll987l ). In the model dis- 
cussed here 7 C is about a few and the growth rate (1221) or 
P3l> can be quite fast. For 7 -> j c = 2, M L = 0.1, 
gives F/lo m 0.1. Although the growth rate can be reduced 
due to a spread in the particle momentum distribution, it 
seems that the instability can develop rapidly. It should be 
pointed out here that radio emission from two-stream in- 
stability is a two-stage process. The instability can generate 
Langmuir waves which can in turn be converted to electro- 
magnetic waves and at the same time facilitate the particle 
precipitation. 

Cyclotron/synchrotron maser 

Since particles with small pitch angles escape through 
the loss cone, particle precipitation creates an inversion 
in the pitch angle distribution of th e upward flow- 
ing particles, driving cyclotron maser (|Melrose fc Dulkl 
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19821: iHewitt. Melrose, Ronnmark 11982) o r sync h rotron 
maser dZheleznvakov fc Suvo'rovf Il972l lYoonl Il990l ; 
iHoshino fc Aronsl Il99lh " This mechanism was discussed 
widely in interpretation of planetary radio bursts (Melrose 
19931 ; lErgun, et af1l2000t ) and solar bursts l|Melrose fc Dulk 
1982 ). Synchrotron maser is applicable in the relativis- 
tic regime p± > 1, while cyclotron maser may be the 
more relevant when such inversion occurs close to the 
star during the precipitation when their gyration is in 
the nonrelativistic regime. The relevant radio wave can 
grow when the frequency satisfies the cyclotron resonance 
condition. For cyclotron resonance at the fundmental 
harmonic uj — Q. e H — k\\V\\ = 0, the growth rate is given by 
jMelrose fc Dulklll982l ) 



(24) 



where Nb is the number density of particles in cyclotron res- 
onance, and i>o is the velocity of particles in cyclotron reso- 
nance. For a typical pulsar with P = 1.5 s and B s = 10 9 T, 
this condition can be satisfied at a radius O.IRlc- It is worth 
noting that the related but different process — cyclotron in- 
stability due to the anomalous doppler effect, which has 
been considered in the literature for pu lsar radio emission 
l|Kazbegi. M achabcli. & Milli kidzel Il99ll ). is not applicable 
here; it requires an energetic beam (7 ~ 10 B — 10 7 ) of parti- 
cles that satisfy the anomalous doppler condition. 



6.3 Propagation 

Propagation of radio emission in the CFLR is constrained by 
induced scattering by electrostatic wav es excited by the ra- 
dio emission itself l|Luo fc Melr ose 2006). Assuming a bright- 
ness temperature T& and a beaming s olid angle Afi, one ma y 
define a radius of an opaque sphere (|Luo fc Melr ose 2006), 



3 x 10 



( 



KJ U.4GHzJ \2s) 



(- 



V10 24 

B x 3/r 



-10/7 



(25) 



V10 9 T, 

such that the wave with a frequency v can only propagate 
outside the sphere 77 > rj cr . Here we assume that the plasma 
density in the propagation path is the local GJ density. 

Although AQTb is not well constrained, the condition 
(|25[) is quite robust. For an observed flux density S° " and 
a pulsar distance Dl, the brightness temperature times the 
beaming solid angle can be estimated from 

2 '1.4GHz\ 2 



AQT b w 2.3 x 10 23 K ( A " 



Dr 



3 kpc 



V10 6 m 

^ f gobs 



) m~ 



(26) 



where As is the distance between the scattering region and 
the emission region. In the case of inward emission due to 
two-stream instability, (|25[l gives a lower limit to the emis- 
sion radius; for the emission to b e potentially visible, the 
emission radius must be 77 > ^r\ cr (|Luo fc Melrosell2006l ). 



6.2 Beaming properties 

Bursts produced from particle precipitation are character- 
ized by sharp spikes. In the case of two-stream instability, 
since emission arising from the instability is confined to a 
thin shell region in the CFLR near the last open field lines, 
one expects the beaming direction of the radio emission to 
be similar to that of reverse emission in the conventional 
polar cap models (|Dvks. et all 120051 ; IZhang. Gil fc Dvksl 
2007). Consider a trapping region has a radial range from 
Rlc(1 — Ag) to Rlc- The inward tangent angles at 77 are 
respectively given by ip2 and ipi. One has a conal structure 
with an angular thickness as Aij) = tp2 — tpi fa (3/4)A£Wi 1//2 
for Ag <C 1. For the emission to be in radio, the emission 
radius has to be 77 ~ 10~ 2 , giving Atp « 0.05 for Ag = 0.1. 
Thus, the model predicts a rather narrow burst profile. Since 
particle precipitation and subsequent burst emission can be 
highly localized and occurs over a large range of altitudes, 
such bursts can appear on a wide range of pulsar phase. 

For cyclotron maser, the beam width is A6 ~ vo /c <C 1 
IIHewitt. Melrose. Ronnmarklll982f ); thus, the maser emis- 
sion can produce narrow spiky profiles. Since cyclotron 
maser emission occurs near particle reflection, the emission 
is beamed at a large angle to the field lines. This beaming 
feature is very different from the usual relativistic beaming 
(I/7 >C 1) along the field lines. Such unsual beaming fea- 
ture can be tested observationally if both bursts and weak 
emission (i.e., the usual pulsed emission) can be detected 
and relative phase between the two components can be de- 
termined. 



7 CONCLUSIONS AND DISCUSSION 

We consider possible existence of transient radiation belts in 
the magnetospheres of pulsars with low LC magnetic fields. 
It is suggested that particle acceleration in latent outer gaps 
can lead to creation of pairs with large pitch angles in the 
CFLRs near the LC where they are trapped due to the 
magnetic mirror effect. In the trapping regions where the 
magnetic fields are weak (~ 10~ 4 T), particles radiate away 
their perpendicular energy over a time much longer than a 
typical pulsar period. Thus, the plasma density can build 
up over a synchrotron cooling time, which is much longer 
than the pulsar period. Similar to the planetary radiation 
belts, pulsar radiation belts can be disrupted due to distur- 
bances of waves propagating into the regions. Catastrophic 
disruption of these regions may lead to intense particle pre- 
cipitation, which in turn generates radio bursts that can be 
seen as radio transients. The main triggering mechanism is 
scattering of trapped particles into a loss cone as a result of 
interactions with waves generated in plasma instabilities in 
the trapping regions or Alfven waves emitted at the surface 
due to low-intensity starquakes or oscillations. The plasma 
that rushes through the background corotating plasma can 
produce bursts of coherent radio emission through stream- 
ing instability or cylcotron/synchrotron maser. Existence of 
such radiation belts, if observationally confirmed, would im- 
ply that some of the bursty phenomena of pulsar radio emis- 
sion may be due to particle precipitation. In particular, we 
suggest that such transient radio emission may be seen as 
RRATs. 

A notable feature of our model, as applied to RRATs, is 
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that in comparison with conventional polar cap models, the 
transient emission mechanism considered here is less con- 
strained by the pair production efficiency; it predicts that 
pulsars below the 'deathline' can still produce radio bursts. 
One of the major problems with polar cap models is the 
low efficiency of pair production. Most RRATs have long 
periods and for pulsars with such long periods, there is in- 
sufficient supply of pairs needed for coherent emission. Al- 
though an outer gap can supply additional pairs t hrough 
accel e ration of charged particles supplied externally ( Chend 
1 19851 ; iRuderman fc Chend Il988l ; ICordes fc Shannon! l200fC 
the pair production efficiency is too low to produce a sub- 
stantial downflow of pairs; the only possible location for re- 
verse emission in the polar regions is near or in the polar 
gap where a downward pair cascade can occur; for example, 
externally-supplied particles flush the polar gap resulting in 
a down ward transient pair cascad e followed by quenching of 
the gap JCordes fc Shannonll2006h . However, such a scenario 
would predict an emission region very close to the star and 
any reverse emission may be eclipsed by either the star or a 
dense plasma near the surface. In our model, since particle 
precipitation occurs in the CFLR, it does not quench the 
polar gap, allowing the possibility of observing both burst 
and normal radio emission. 

In application to magnetars with active high-energy 
transient emission, one expects low-intensity quakes to 
occur much more frequen tly than in normal pulsars 
l|Thompson fc Duncanlll99"5l ) . Such frequent disruption may 
prevent a radiation belt forming in a magnetar magneto- 
sphere, even though the synchrotron decay time is relatively 
long near the LC. This scenario seems to be supported by 
lack of any detection of RRAT-like radio bursts from the 
known magnetars. Although transient radio emission was 
recently detected from the magnetar XTE J1810— 197, po- 
larization study suggests t hat the emission geo metry resem- 
bles that of young pulsars (|Camilo et al.ll2007h . 

In principle, our model may also apply to old, long- 
period pulsars. The loading time tl can be significantly 
longer than the synchrotron decay time. The plasma trap- 
ping can be due to in situ ionization of the accreted neu- 
tral grains in the CFLR. Since the synchrotron decay time 
is long, for example, t s « 16 yr for PSR J2144— 3933 with 
P = 8.5 s and B s — 4 x 10 8 T, a nearby long-period pulsar, 
accumulation of plasma can occur over a long time (provided 
that the reoccurrence time of the triggering mechanism is 
also long). Since old pulsars are not energetic, the predicted 
radio luminosity is low. As it takes a long time to replenish 
the trapped plasma, such transient events can be rather in- 
frequent but may still be detectable provided these pulsars 
are located relatively nearby. 
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APPENDIX A: BOUNCING TIME 

From the first adiabatic invariant p±/B = const and the 
energy conservation p 2 — const, one has 

2 B( V ) 2 



B(r)i) 



■P\\(.V), 



(Al) 



where the relevant quantities are functions of r\. Assuming 
sin a; = p±(jji)/p, the time to travel from the injection ra- 
dius rji to the reflection radius rj r and back to r\i can be 
written in terms of integration along the particle's path ds: 

2 f ds 



At 



Rlc 



[1 - B( V )/B( Vr )}^ 2 

1/2 



4-37; 



dr) 



V J ll-(Vr/v) 3 ] 1/2 ' 



(A2) 



where we assume a particle follows the last closed field 
line of a dipole field and r\ T < r\i satisfies the condition 
1 — [B(rj r )/ B(r]i)] sin 2 on = 0. The integration yields ap- 
proximately At « rjiP/TT. One may assume without loss of 
generality r)i = 1 to show that the bounce time (back and 
forth between two mirror points) is — 2At ~ 2P/n. The 
bounce time is not sensitive to the pitch angle at the injec- 
tion provided that the condition on B\ is satisfied; this 
condition requires the reflection radius to be larger than the 
stellar radius, i) r ^ Rq/Rlc- 



(B3) 



V(s,p,k) 



(±v ± [e»J,- l {z)+e- it '*J a+ i(z)] , 
-~i5 e v ± [e**J-i(z) - e- is °*J s+1 (z)] , 
v\\Js(z)), (B4) 



with S e = 1 for positrons and <5 e = — 1 for electrons, R is 
the ratio of the electric to total wave energy density, z — 
k±v±y /H e and k = (k± cos0, k± sin<^>, fey). In the following 
we consider the small 2«1 limit. The Bessel function has 
an approximation J s ~ (z/2) s /V(s + 1) for s ^ 0. This gives 
Jo = — Ji an d J_i = J\. 

The polarization of a low-frequency Alfven wave can be 
written as 

e A w (cos 0, sin 0, e A \\), (B5) 

where e A \\ ~ (u> lu> v ) 2 ( k_i_/k\\)n (lArons fc Barnard! 1 19861 : 
iMelrose fc Gedalini ri999), and n is the refractive index. It 
can easily be verified that e A ■ V = e^iiUM for s = 0. Assum- 
ing R = 1/2 and sin# = k±/k, one obtains 



w(0,p,k) 



eoTiu) 



(B6) 



The diffusion coefficient is 



where uj « knc, Ub = SB A /2fio is the magnetic energy den- 
sity of the wave concerned, and U p = m e c 2 (^NlMl is the 
density of the plasma kinetic energy. It should be empha- 
sized here that the disffusion in pitch angle is caused by 
change in pu not p±. 



APPENDIX B: PITCH-ANGLE SCATTERING 



We outline quasilinear theory for pitch-angle scattering. Let 
f(p, a) be a gyrophase-averaged distribution in the particle 
momentum space. In the quasilinear diffusion formalism one 
has 



dt 



1 



d 



sin a da 



p 2 dp 



Dpa da +Dpp d P ' f 



(Bl) 



where the diffusion coefficients are written in terms of the 
scattering pro bability w(s,p , k) and the wave occupation 
number N(k) l|Melrosdll98oh : 











Dpp _ 
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dk 

j^-^w{s,k,p)N(k) 



(Aa) 2 
AaAp 
L (Ap) 2 



(B2) 



with D ap = D pa , Aa = h(uicosa — k«v)/(pvsiaa), and 
Ap = Tkv/v. The scattering probability is 

w(s,p,k) = 2 ^ h ^ \e ■ V(s,p,k)\ 2 



